clear;clc
x=0:0.1:1;
fx=x.^6/sqrt(6+x)+exp(x);
trapz(x,fx)
x=linspace(0,5,200);
y=x.^6/sqrt(6+x)+exp(x);
plot(x,y)
p=polyfit(x,y,5);
y1=polyval(p,x);
hold on
plot(x,y1)
legend({'actual','fitiing'})